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One of the non-Abelian features of parton energy loss is the ratio AE g / AE q = 9/4 between 
gluon and quark jets. Since jet production rate is dominated by quark jets at high xt = 2pr/y/s 
and by gluon jets at low xt, high pr hadron suppression in high-energy heavy-ion collisions should 
reflect such a non-Abelian feature. Within a leading order perturbative QCD parton model that 
^j- ' incorporates transverse expansion and Woods-Saxon nuclear distribution, the energy dependence of 

large pr ~ 5 — 20 GeV/c hadron suppression is found to be sensitive to the non-Abelian feasture of 
parton energy loss and could be tested by data from low energy runs at RHIC or data from LHC. 

(N 
o 

One of the ultimate goal of the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory is 
to produce the quark-gluon plasma (QGP) by smashing two gold nuclei at the speed of light. The discovery of 

0, |j| in central Au+Au collisions together with the observation of parton recombination 
and the early thermalization of the dense matter [H [3JTl|ll Il7l 111 liH I20L I2I 15^ has 
provided clear evidence for the formation of strongly interacting partonic matter~|23l l24j . The observed jet quenching 
effect manifests itself in several aspects of high px hadron spectra, which include suppression of inclusive spectra in 
■^-j- [ Au+Au relative to pp collisions pl@]j disappearance of back-to-back correlation Q and the azimuthal anisotropy in 
non-central Au+Au collisions Q. The absence of these jet quenching phenomena in d + Au collisions |25ll2rl I27I |2S| 
shows that they are due to final state interactions with the produced strongly interacting matter. Detailed analyses 
indicate that parton energy loss is the source for the observed jet quenching [2^, |3(], |32, H3, l34j . The initial 
gluon density, which the parton energy loss is proportional to, has been extracted from RHIC data of central Au+Au 
collisions at y/s = 200 AGeV and is about 30 times higher than that in a cold nucleus |35|, |36( . 

The radiative parton energy loss incoporated in previous studies within a leading order (LO) perturbative QCD 
(pQCD) model [Hi El M, El has two basic non-Abelian features. One of them is the quadratic dependence 
on the total distance traversed by the propagating parton due to the non-Abelian Landau-Pomeranchuk-Midal (LPM) 
interference effect in gluon bremsstrahlung induced by mulitple scatterings in a static medium 0, 0> El lil I45L Ha . 
[H Eflj . The second feature of the parton energy loss is its dependence on the color representation of the propagating 
parton. Therefore, energy loss for a gluon is 9/4 times larger than a quark. Previous works have investigated the 
consequences of the second non-Abelian feature in the flavor dependence of the high-pT hadron suppression |38| . In 
' this paper we study the effect of the non-Abelian parton energy loss on the energy dependence of the inclusive hadron 
spectra suppression. We exploit the well-known feature of the initial parton distributions in nucleons (or nuclei) 
that quarks dominate at large fractional momentum (x) while gluons dominate at small x. Jet or large pt hadron 
production as a result of hard scatterings of initial partons will be dominated by quarks at large xt = 2pt/\/s and 
by gluons at small xt- Since gluons lose 9/4 times as more energy as quarks, the energy dependence of the large (and 
fixed) pt hadron spectra suppression due to parton energy loss should reflect the transition from quark-dominated 
jet production at low energy to gluon-dominated jet production at high energy. Such a unique energy dependence of 
the high-pT hadron suppression can be tested by combining y/s = 200 AGeV data with lower energy data or future 
data from LHC experiments. 

We will work within a LO pQCD parton model incorporating the non-Abelian QCD parton energy loss in high- 
energy heavy-ion collisions. We will study the energy dependence of the high-py hadron suppression and compare the 
effect of QCD energy loss with that of a non-QCD one where gluons and quarks are chosen to have the same amount 
of energy loss. In both cases, we will assume that parton energy loss is proportional to the initial gluon density of the 
system which in turn is assumed to be proportional to the measured total charge hadron multiplicity in the central 
rapidity region. 

In comparison to previous studies within the LO pQCD parton model that employed the hard-sphere model of 
nuclear distribution and assumed only longitudinal expansion, we will use more realistic Woods-Saxon nuclear distri- 
bution and in addition include the tranverse expansion of the dense medium. 

In a LO pQCD model 37J , the inclusive invariant differential cross section for high-py hadrons in A + B collisions 
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where a(ab — > cd) are elementary parton scattering cross sections. The factor K w 1.0 — 2.0 is used to ac- 
count for higher order QCD corrections and is set to be the same for both p + p and A + B collisions at the 
same energy. The hadron is assumed to have the same rapidity as the parton, i.e. y — y c , and its frac- 
tional momentum is defined by z c — pt/ptc- The parton distributions per nucleoli f a /A{ x a,Q 2 ,f) inside the 
nucleus can be factorized into the parton distributions in a free nucleon given by the CTEQ parameterization 
|49l l50| and the impact-parameter dependent nuclear modification factor given by the new HIJING parameteri- 
zation: f a/A (xa,Q 2 ,r) = R^(x,Q 2 )[(Z/A)f a/p ( Xl Q 2 ) + (1 - Z/A)f a/n (x,Q 2 )] with R^x^Q 2 ) given by Eqs. (8) and 
(9) of Ref. [51j. We assume that the initial transverse momentum distribution gA(kx,Q 2 ,b) has a Gaussian form 
|3jj|52j with a width that includes both an intrinsic hr in a nucleon and the nuclear broadening due to initial multiple 
scattering in a nucleus: gA(kr,Q 2 ,b) = e~ kr/ A / (7r(fcf,) a)- The impact-parameter dependent broadened variance 
is given by (k^AiQ 2 ) = {k 2 -) n{Q 2 ) + S 2 (Q 2 )[vA{b) — 1], where the number of scatterings v A (b) the projectile suffers 
inside the the nucleus is VA{b) = er/vjv^A^) with the nuclear thickness function t A (b) defined as follows, and the scale- 
dependent S 2 (Q 2 ) is chosen as <5 2 (Q 2 ) = 0.225 ln 2 (Q/GeV)/[l + ln(Q/GeV)] GeV 2 /c 2 . The average initial intrinsic 
transverse momentum in nucleon-nucleon collision is (fcy)jv(Q 2 ) = 1.2 + 0.2Q 2 a s (Q 2 ). The scale which characterizes 
the partonic process is chosen to be Q — px, where pt is the transverse momentum of the final-state partons in a 
partonic scattering. Detailed description of this model and systematic comparisons with experimental data can be 
found in Ref. |3j|. In this paper we use the Woods-Saxon nuclear distribution Fws(r) — Na/[1 + exp((r — Ra)/o)\ 
to replace the simplified hard-sphere one as used in previous papers. Here Ra is the radius of the nucleus and given 
by R A = 1.12A 1 / 3 - 0.86A- 10 / 3 , a = 0.54 fm is a radius parameter and Na is the normalization constant. It can 
be further written as a function of the coordinate component z along the beam direction of the nucleus and b that is 
perpendicular to it by r — V ' z 2 + b 2 . The nuclear thickness function tA(b) is then tA(b) = f_ dzFws(z, b) with the 
normalization condition J d 2 btA(b) — A. 

The parton energy loss is encoded in an effective modified fragmentation function [H^L IH3 | 



D h/c (z c ,Q 2 ,AE c ) - (1 



^M^ + i^fD^Q 2 ) 



D° h/c (zcQ 2 ) 



(2) 



This effective form is a goo d approximation to the actual calculated medium modification in the multiple parton 



scattering formalism 
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given that the actual energy loss should be about 1.6 times of the input value in the 



above formula. Here z' c — Pt/{ptc — AE C ), z' g = (AL / X)pr / AE C are the rescaled momentum fractions and AE C is 
the total energy loss during an average number of inelastic scatterings (AL/ X). The fragmentation functions in free 
space D^ c (z c , Q 2 ) are given by the BBK parameterization |57j . 

In contrast to previous calculations where only longitudinal expansion was considered, we incorporate in this paper 
both longitudinal and transverse expansion of the medium in the calculation of parton energy loss. To simplify the 
calculation, we use hard-sphere nuclear distribution again. Let us denote the gluon number N g and assume that 
N g is a slowly varying function of rapidity y and proper time r at central rapidity region y — 0, then we have 
d 2 N g /drdy = 0. Noting that dN g /dy = pdV/dy and dV/dy 
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The radius has the form Rt{t) = Ra + ( r — T o) c2 s where c s is the speed of sound in the medium given by 
(1/3 for ideal gas). The solution of the above equation is then 

rp [R A + {t- to)c 2 ] 2 = T Q p R 2 A . 



(3) 

(4) 
dP/de 
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The expansion is characterized by the gluon density p g (r, r) whose initial distribution is proportional to the transverse 
profile of participant nucleons. We can write the total energy loss for a parton traversing the medium as 

dE [ T ™* r [R min + (r - t )c 2 s ] 2 - t RI 



A£?(6,r,0) « C-) ld / dr " SJ "~^P g {r, b, r + nr) , (6) 

where R m m = Min(_R J 4, Rb), n is the direction where a parton is propagating. The upper limit T max = Min(AL,r/) 
is the longest time for the parton to propagate in the dense medium, where Tf is the lifetime of the dense matter 
before breakup. AL(b, r, <fi) is the distance the parton, produced at r, travels along n at the azimuthal angle <fi relative 
to the reaction plane in a collision with impact parameter b. Since the formation time of a hadron fragmented from 
a parton is proportional to the energy of the parton, very high energetic partons generally hadronize after the dense 
medium breaks up, or they hadronize outside the medium. In this case we have r max = AL. (dE/dL)id is the average 
parton energy loss over a distance Ra in a 1-dimensional expanding medium with an initial uniform gluon density po- 
The gluon density p g in the longitudinally and transversely expanding medium is then given by 

o„(r,6,r + nr) = - T ° P ° 



[Rmin + (r - Tq)cI] 2c AB Rmin 

(7) 
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The average number of scatterings along the path of parton propagation is 

(AL//) = [ drap g (T,b,r + nr) . (8) 

J T(, 

The energy loss function can be parameterized as 

dE (E/po - 1.6) 1 - 2 

{ dL )ld = 60 7.5 + E/p (9) 

according to a study of parton energy loss 01 that include both included bremstrahlung and thermal aborption of 
gluons. For y / s=200 AGeV, we find following set of parameters en = 1.2 GeV, po — 1-6 GeV and Zo = 0.2 fm (Z 
appears in the formula of (AL/A)) can fit the data. In Ref. [35l Isfll]. these parameters are set to slightly different 
values e = 1-07 GeV, p = 1.5 GeV and Z = 0.3 fm. The value of (|f ) ld with e = 1.2 GeV, = 1.6 GeV 
used in this paper is almost the same in the energy range E = 5 — 20 GeV as with previous values en = 1-07 GeV, 
= 1-5 GeV |Hi|3g|. For example, at E = 5 and 20 GeV we have ) u {e Q = 1.2, ^ = 1-6) = 0.19 and 1.05, while 
(4§)i<i(eo — 1-07, po — 1-5) = 0.19 and 0.99. Another modified parameter Iq in this paper is inversely proportional to 
the average number of scatterings undergone by the propagating energetic parton. The smaller value Zo = 0.2 fm than 
previously used one Zo = 0.3 fm means that the average number of scatterings is tuned larger to make more energy 
loss by compensating the effect caused by transverse expansion which makes the medium more rapidly diluted. Note 
that the parameter eo is proportional and Zo is inversely proportional to the gluon or multiplicity density per rapidity. 
The energy loss in a corresponding static medium is found to be 14 GeV/fm, which is about 30 times as high as in a 
cold nuclei |35| . 

The jet quenching effect can be shown by the nuclear modification factor defined as [58| 

R AB = d<J AB/ d yd 2 PT , 1Q , 

(^binary 

where -/Vbmary is the average number of geometrical binary collisions at a given range of impact parameters 

(AW y ) - J d 2 hd 2 vt A {r)t B {\h~r\) . (11) 

If there is no energy loss, the cross section for nucleus-nucleus collisions is a simple sum of that for elementary binary 
nucleon-nucleon collisions, so the nuclear modification factor Rab is one. Hadron suppression due to parton energy 
loss leads to Rab < 1- 

As we mentioned earlier that we use the Woods-Saxon nuclear distribution in the parton model calculation. The 
numerical difficulty with the Woods-Saxon distribution is that one cannot simply put the analytical formula into 
the program becasue that would substantially slow the speed of the calculation and make the numerical calculation 
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practically impossible. One trick to overcome this problem is to calculate the distribution before hand and then store 
the results in tables whose entries can be called in the run time of the program. The calculated Rab results for Au+Au 
collisions are shown in Figs. 1 and 2 with Fig. 1 for ir° and Fig. 2 for charged hadrons. The results for three collision 
energies -^=62.4, 200 and 5500 AGeV are given. The parameters eo and Ao at these energies are set to approriate 
values based on the ratios of charged particle (or gluon) multiplicity density [^1 EHI (dN c h jdy) 5500 /(dN c h/dy) 200 , 
(dN c h/dy)s2A/ {dN c h/dy)2oo and their values eo = 1.2 and Ao = 0.2 at 200 AGeV. In the figures one can see different 
transverse momentum behaviors of the nuclear modification factor at these energies. Similar behaviors have been 
seen in recent studies 39, 59, 60]. The nuclear modification factor decreases with px at 62.4 AGeV, while it slightly 
increases with px at 200 AGeV. So the nuclear modification factors for neutral pions and charged hadrons at 62.4 
AGeV intersect at about px = 11 and 10 GeV respectively with those of 200 AGeV in the QCD case, where the 
energy loss parameters for gluons and quarks satisfy g/q — 9/4 (we will explain this point later). The same feature 
also occurs in Ref. |39j| where the hard sphere distribution and only the longitudinal expansion are used. In the 
intermediate px region, one expects the jet fragmentation process to be modified by other non-perturbative processes 
such as parton recombination or coalescence jSj, 0, |n| . The observed flavor dependence of the hadron suppression 
and of the azimuthal anisotropy clearly points to the effect of parton recombination that enhances both baryon and 
kaon spectra in the presence of dense medium. To include this effect in the current parton model, we have added a 
soft component to kaon and baryon fragmentation function that is proportional to the pion fragmentation function 
with a weight ~ (N\yi n (b,r))/[1 + exp(2px c — 15)] where px c is the transverse momemtum for parton c (actually we 
have also found that the similar effect can be also achieved by using a function of the hadron transverse momemtum 
Pt- (Nbm(b, r))/[l + exp(pT — 5)].) The functional form and parameters are adjusted so that (K +p)/n « 2 at px ~ 3 
GeV/c in the most central Au+Au collisions at y/s — 200 AGeV and approaches its p + p value at px > 5 GeV/c. 
This gives rise to the splitting of the suppression factor for charged hadrons and ir° in the calculation. 

To study the sensitivity of hadron spectra suppression to the non-Abelian parton energy loss, we compare the 
results with two different cases at each energy: one for the QCD case where the energy loss for a gluon is 9/4 times as 
large as for a quark, i.e. AE g / AE q = 9/4, the other is for a non-QCD case where the energy loss is chosen to be the 
same for both gluons and quarks. Similarly, the average number of inelastic scatterings obeys {^ L )g/{^ ± )q = 9/4 
in the QCD case. For the non-QCD case we are considering, the above ratio is set to one. In Figs. 1 and 2 one 
can see that the difference between the QCD and non-QCD cases are more significant for higher collision energies. 
This fact manifests itself at 200 and 5500 AGeV, where the nuclear modification factors Rab are much lower for the 
QCD energy loss pattern than for the non-QCD one. As shown in the figures, the suppression at 62.4 AGeV is not 
sensitive to gluon energy loss, but only to quark energy loss because of the dominance of quark jets at large pr- At 
200 AGeV, however, the suppression is sensitive both to quark and gluon energy loss. At LHC energy, the results 
are only sensitive to gluon energy loss in the px range we calculated. Such energy dependence pattern is a direct 
consequence of the non-Abelian feature of the energy loss. 

In order to demonstrate the colliding energy dependence of the nuclear modification factor and illustrate the 
difference between QCD and non-QCD energy loss, we compute the Raa for neutral pions at fixed pr = 6 GeV in 
Au+Au collisions as a function of y/s from 20 AGeV to 5500 AGeV. Shown in Fig. 3 are the calculated results with 
both the QCD energy loss and a non-QCD case where the energy loss is set to be identical for quarks and gluons. 
Two parameters eo and lo which are relevant to the energy loss are determined according to the gluon number or 
the charged particle multiplicity per rapidity |5lL l6l| . One can see that due to the dominant gluon bremstrahlung or 
gluon energy loss at high energy the Raa for the QCD case is more suppressed than the non-QCD case where the 
gluon energy loss is assumed to take an equal role as the quark one. In the calculation, we have assumed that the 
lifetime of the dense matter is equal or longer than the parton propagation time which is essentially determined by 
the system size. This might not be the case for heavy-ion collisions at lower energies, in particular at around y/s = 20 
AGeV. If one takes short lifetime, the suppression factor Raa is much larger than 1 due to strong Cronin effect [39]] . 
The dashed box around y/s = 20 AGeV in the Fig. 3 assumes a lifetime r/ =0 — 2 fm/c and thus provides an estimate 
of the uncerntaintty due to lifetime of the dense matter. Since finite lifetime reduce the effect of full parton energy 
loss, the difference between QCD and non-QCD energy loss effect in Raa should be smaller. The difference in Fig. 3 
is therefore the upper limit. 

Another interesting feature with the energy dependence of Raa is the change of slope around y/s — 1300 GeV. The 
rapid decrease of Raa at y/s — 20 — 1300 GeV is mainly due to increase of parton energy loss due to increased initial 
gluon density and also the change of px slope of jet production cross section with y/s. As the energy loss increases, 
more jets produced inside the overlapped region are completely suppressed. Only those that are produced within an 
outlayer in the overlapped region will survive. This will be like surface emission with a finite depth. The suppression 
factor Raa will then be determined by the width of the outlayer which is just the averaged mean-free-path (A). As a 
consequence, Raa will then have much weaker y/s dependence. 

In summary, nuclear modification factors in Au+Au collisions at -^5=62. 4, 200 and 5500 AGeV are calculated in 
a LO perturbative QCD model with medium induced parton energy loss. The previous calculations based on the 
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FIG. 1: Nuclear modification factor RauAu for neutral pions at y^s =62.4, 200 and 5500 AGeV. We choose the values of 
corresponding parameters at 62.4 and 5500 AGeV based on their values at 200 AGeV and the ratio (dN c h/ dy)^2 / {dN c h / dy)2oo 
and (dN ch /dy) 550 o/(dN ch /dy) 2 oo- 



hard-sphere distribution of nucleus and the longitudinal expansion of the dense medium are improved in terms of a 
more realistic Woods-Saxon distribution and both longitudinal and transverse expansion. The comparison of nuclear 
modification factors for energy loss patterns in QCD and non-QCD cases shows sizable difference at higher colliding 
energies and thus can be tested by the energy dependence of the hadron suppression factor Raa in the range between 
\fs = 20 — 1000 GeV. We also found a weaker energy dependence above y/s = 1000 GeV due to surface emission with 
finite depth. 
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FIG. 3: Nuclear modification factor RauAu for neutral pions as function of collision energy at fixed pr — 6 GeV in most central 
collisions (with centrality 10%). Here we compare the QCD energy loss and a non-QCD one where the energy loss is identical 
for quarks and gluons. 
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